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Abstract 

In this paper, we present a model for determination of a weakly time dependent effective recombination coefficient for the perturbed 
terrestrial ionospheric D-region plasma. We study consequences of a class Ml. X-ray solar flare, recorded by GOES- 15 satellite on 
February 18, 2011 between 14:00 UT and 14:15 UT, by analyzing the amplitude and phase real time variations of very low frequency 
(VLF) radio waves emitted by transmitter DHO (located in Germany) at frequency 23.4 kHz and recorded by the AWESOME 
receiver in Belgrade (Serbia). Our analysis is limited to ionospheric perturbations localized at altitudes around 70 km where the 
dominant electron gain and electron loss processes are the photo-ionization and recombination respectively. 

^ Keywords: electron concentration, photo-ionization, electron recombination, solar flare, ionosphere 
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1. Introduction 

\ The ionosphere, having characteristics of plasma, is very sen- 
sitive to electromagnetic disturbances whose intensity and num- 
ber vary with solar activity. These disturbances cause numerous 
complicated physical, chemical and dynamical phenomena in 
the lower ionosphere and may directly affect human activities, 
especially in the telecommunications industry. Besides a pure 
Scientific interest to study the influence of solar activity on the 
terrestrial atmosphere, the understanding and predicting the re- 
,sulting turbulent regions of the ionosphere has important appli- 
cations for radio communications, military operations in remote 
locations, planned networks of mobile communications satel- 
lites, high-precision applications of global navigation satellite 
Isystems, etc. 

' The VLF radio signal propagation properties are determined 
by the wave attenuation and reflection which depend on the 
electron concentration and collision frequency. The elec- 
■ tron concentration calculations in empirical and semi-empirical 
models are based mainly on measurements using radio-wave 
propagation methods and sounding rockets. 

Due to a varying complex structure of the atmosphere and be- 
cause of numerous influences coming from the Earth and outer 
space, experimental simulations of processes in the ionosphere 
are very difficult. For this reason, theoretical models are of a ba- 
sic importance in providing us with information on ionospheric 
parameters, with the effective recombination coefficient a^s be- 
ing one of them. This coefficient plays an important role in the 



ionization balance and in kinetics of underlying chemical reac- 
tions and from its values one can estimate the structure of the 
lower ionosphere. 

In the Earth ionosphere, the electron concentration depends 
on photo-ionization and recombination processes. At heights 
around 70 km, the dominant recombination processes are the 
electron-ion, ion-ion and three-body recombination character- 
ized by a common effective recombination coefficient a^s fl] . 

In the case of unperturbed, quasi-stationary ionosphere, the 
electron gain and electron loss processes are in equilibrium, i.e. 
mutually balanced. 

During a solar flare, the electron concentration is raised due 
to increased photo-ionization and this time period can be named 
as photo-ionization regime. After a flare maximum when the 
rate of recombination becomes larger than the photo-ionization 
rate, the electron concentration starts to decrease and the iono- 
sphere enters a recombination regime. 

In this paper, we calculate the time dependence of effec- 
tive recombination coefficient (Xeif for different altitudes in the 
ionospheric D-region for the tail of the recombination regime. 
The considered class Ml.O X-ray solar flare was recorded by 
GOES-15 satelHte on February 18, 2011 between 14:00 and 
14: 15 UT. Our analysis uses a discrete set of data on amplitude 
and phase variations of the 23.4 kHz VLF signal emitted by 
the transmitter DHO in Germany and registered by the AWE- 
someQ receiver system H in Belgrade Serbia with the sam- 
pling period of 0.02 s. These experimental data were treated 
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by the LWPC (Long-Wave Propagation Capability) computer 
program for modeling the ionosphere which yielded dis- 
crete sets of data on the wave reflection heights H'{t) and sharp- 
ness Pit) the basic parameters of Wait's model of the ionosphere 
101 needed to compute the distribution profiles of electron con- 
centration N(t,h) and the effective recombination coefficient 
o^esit, h) where h is the location altitude. 

2. Basic theory 

The electron concentration N{t,h) at altitude h in the D- 
region is time dependent and is taken to obey the following 
differential equation [5]: 

dN(t, h) 



dt 



= q{t, h) - ffetf (f, h)N-^{t, h), 



(1) 



where the concentrations of negative ions and their time deriva- 
tives are comparatively small at altitudes of about 70 km L6J so 
that their contributions are not taken into account. 

According to the Chapman's theory ff\ the rate of electron 
production q{t, h) can be expressed as: 



(2) 



where I(t, h) is the X-ray radiation flux intensity and Q\{t, h) is 
a proportionality coefficient whose explicit expression turns out 
not to be required in our computations. 

The X-ray radiation flux intensity I(t, h) at some level h is 
related to its value lao (measured by the satellite GOES-15) 
through the following relation: 

I{t,h)^Q2(t,h)Ut), (3) 

where Q2{t,h), like Q\{t,h) in Eq. (|2]l, is a proportionality co- 
efficient whose explicit expression will not be needed in our 
further treatment. 

From Eqs ([TJ-(l3]l we can write: 

dN(t, h) 



dt 



K(t, h)IUt) - tteif (f, h)N^(t, h). 



with the coefficient K{t, h) given by: 
K{t,h) = Qi(t,h)Q2(t,h). 



(4) 



(5) 



After a sufficient time interval following the intensity maxi- 
mum (the shaded domain in Fig.[T]), the coefficients K{t, h) and 
Oesit, h) become weakly time dependent and we will consider 
them approximately constant within some finite time interval 
At: 

K{t - At, h) K Kit, h) = Kit, h), 

Uesit — At, h) K Ogfiit, h) = aeffit, h). 

Eq. Q can now be applied to the interval end points fi = 
t — At and t2 = t which yields the following set of two algebraic 
equations for unknown quantities Kit, h) and Oesit, h): 

dN — 

^ Kit,h)IUt-At) 

t-Ai,h 

-aeffit, h)N^it - At, h) 
= Kit, h)IUt) - aeffit, h)N^it, h) 



dt 



(6) 



dt 



which finally gives the expression for numerical computation of 
the weakly time dependent effective recombination coefficient 

aeffit, h): 



aesit, h) 



looit - At) dN 
Jiit, At, h) dt 



Ut) dN 



Jlit, At, h) dt 



t-At,h 



(7) 



where: 



Jlit, At, h) = N^it - At, h)I^it) - N^it, h)IUt - Af). 



(8) 



All quantities on the right-hand-side in Eq. O have values 
that can easily be obtained from both the recorded signal prop- 
erties and Wait's model of the ionosphere 101. In computations 
that follow, we take Af = Is. 

3. Experimental data and calculation procedure 

Within the time interval 14:00 UT - 14:15 UT, February 18, 
2011, a solar flux increase was registered in the range 0.1- 
0.8 nm by the GOES-15 satellite whose maximal value cor- 
responded to the level typical of a class Ml.O solar X flare 
(Fig. [TJ. This rise on the radiation flux altered the rate of 
photo-ionization and electron recombination process and, con- 
sequently, it changed the electron concentration. As a result, 
the recorded VLF wave exhibits a time varying amplitude and 
phase as shown in Fig.[T] 

To obtain the effective recombination coefficient aeffit,h) 
profiles given by Eqs. O-® we need the explicit expression 
for the time dependence of the electron concentration profile 

and its time derivative dN/dt as well as the incoming solar 
radiation flux as recorded by the GOES-15 satellite. 

The electron concentration A^(f, h) is computed by the follow- 
ing relation |4]: 



Nit, h) 



1.43 • ioi^^e-'8W«'(')eW)-o.i5);i 



(9) 



where H'it) and (Ut) are defined as the wave reflection height 
and sharpness respectively. Quantities H'it) and Pit) are deter- 
mined by the LWPC computer program |j3j] using the recorded 
discrete time-set of values for the amplitude and phase of the 
considered VLF 23.4 kHz VLF signal (iS §1). Time depen- 
dencies of these two parameters are shown in Fig. |2] where it 
can be seen that the VLF wave reflection height H'it) initially 
decreases in time and reaches a minimum at about two min- 
utes after the solar flux maximum (Fig.|2]l. The sharpness Pit), 
characterizing the gradient of the electron concentration with 
height, has a reversed time distribution, i.e. it first increases 
and attains a maximum at the same time when H'it) has its 
minimum. 

As the recorded quantities are discrete in time, the input ra- 
diation flux intensity /oo(f) and the derived values for Pit) and 
H'it) are not smooth functions of time f and we fitted them by 
a continuous model function implemented in the Origin-8 pro- 
gram, also found in the NAG (Numerical Algorithms Group) 
library: 



t.h 



I + e "2 



1 

^ 1 + e "3 / 



(10) 



where x is time t in seconds starting from 14:04:00 UT and y 
stands for Ioc{t), /3(t), and H'(t). The numerical fitting constants 
yo. A, Xc, wi, W2 and W3 related to the considered physical quan- 
tities y in Eq. ( fTOl l were obtained as the program outputs and are 
presented in Table 1 . 

Thus derived continuous time dependent function N{t,h) 
yields a smooth time derivative dN/dt and the effective recom- 
bination coefficient a^sit, h) according to Eqs. ©-(1911. 

4. Results and discussion 



computer program. In the case of this flare, the effective recom- 
bination coefficient has values within the order of magnitude 

10-12 m^s-i to 10-" m^s-i. 
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The electron concentrations at various heights and their time 
derivatives obtained from Eq. (|9]l are shown in Figs|3]and2]re- 
spectively. We can see that the considered solar flare causes 
larger increases in electron concentration at higher altitudes. 
Comparing Fig. [1] and Fig. |3] it can be noticed that the time 
distribution of the electron concentration follows the time vari- 
ation pattern of the registered solar flux on GOES-15 satellite. 
Maximum of solar flux occurs at 11: 24 UT which is about two 
minutes before the time when the corresponding maximum ap- 
pears in the electron concentration curve. At considered alti- 
tudes of around 70 km, the dominant electron gain and electron 
loss processes are the photo-ionization and recombination re- 
spectively. Within the time period before the electron concen- 
tration maximum, the rate of the photo-ionization is larger then 
that of recombination processes. Contrary, the domination of 
recombination processes causes a decreasing electron concen- 
tration as seen in the recombination regime in Fig. |3] 

The rate of the electron concentration change is more pro- 
nounced at higher altitudes. It is positive when photo-ionization 
processes dominate the recombination processes. Contrary, the 
domination of the recombination produces a time decrease of 
the electron concentration and, consequently, its negative time 
derivative. 

In Fig.[Tl the shaded time domain indicates a relatively slow 
time variation of indicated quantities that allows us to use the 
procedure for calculation effective recombination coefficient 
a'eff(f, h) according to the described basic theory. Time varia- 
tions of aeff(f, h) are presented in Fig. |5]for different altitudes. 
Unlike the electron concentration, the effective recombination 
coefficient has larger values at smaller altitudes. 

The increase of coefficient Oefiit, h) is in agreement with ^ 
where a^sit, h) has larger values for smaller maxima of flare 
radiation intensity. 



5. Conclusion 

In this work, we present a method to calculate the effec- 
tive recombination coefficient at altitudes around 70 km for the 
tail of the recombination regime. It is shown that the effective 
recombination coefficient depends on the solar flux intensity, 
electron concentration and its time derivative. We present a 
method of determining the electron concentration by compar- 
ing the recorded VLF signal amplitude and phase changes with 
the corresponding values obtained numerically using the LWPC 
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Figure 1: Solar flux registered by GOES-15 satellite and perturbed amplitude 
and phase of signal emitted from transmitter DHO (Germany) and recorded on 
AWESOME receiver in Belgrade (Serbia) during observed flares. 
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Figure 2: The reflection heiglit H' and sliarpness fi obtained by comparison 
LWPC simulation and recorded amplitude and phase values. 




Figure 3: The time distribution of electron concentration during photo- 
ionization and recombination regimes. 
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Figure 5: Effective electron recombination coefficient in ionospheric D-region 
for the tail of the relaxation regime following solar flare on February 18, 201 1. 
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